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ABSTRACT 

High resolution X-ray observations indicate that the entropy profiles in the central regions of some 
massive cooling flow clusters are well approximated by powerlaws. McCarthy and coworkers recently 
accounted for this trend with an analytic model that includes the detailed effects of radiative cooling. 
Interestingly, these authors found that cooling (and subsequent inflow of the gas) naturally establishes 
approximate steady-state powerlaw entropy profiles in the cores of clusters. However, the origin of 
this behavior and its dependence on initial conditions have yet to be elucidated. In the present 
study, we explain this trend in the context of the self-similar cooling wave model developed previously 
by Bertschinger (1989). It is shown that the logarithmic slope of the entropy profile in the cores 
of relaxed cooling flow clusters is given by a simple analytic function that depends only on the 
logarithmic slopes of the local gravitational potential and the cooling function. We discuss a number 
of potentially interesting uses of the above result, including: (1) a way of measuring the shapes of 
gravitational potentials of cooling flow clusters (which may, for example, be compared against the 
standard hydrostatic equilibrium method); (2) a simple method for constructing realistic analytic 
cluster models that include the effects of radiative cooling; and (3) a test of the reliability of cooling 
routines implemented in analytic models and hydrodynamic simulations. 

Subject headings: cosmology: theory — galaxies: clusters: general — X-rays: galaxies: clusters 
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1. INTRODUCTION 

The presence of strong positive temperature gradients 
(e.g., Allen et al. 2001; De Grandi & Molendi 2002; Pif- 
faretti et al. 2004; Vikhlinin et al. 2004) and the lack 
of an obvious entropy floor (e.g., Pratt & Arnaud 2002; 
Mushotzky et al. 2003; Piffaretti et al. 2004) in the cores 
of some clusters have stimulated interest in theoretical 
models of the intracluster medium (ICM) that include 
the effects of radiative cooling. Such models have been 
shown to be in broad agreement with the current suite of 
X-ray observations of clusters, particularly if some form 
of feedback/heating is also incorporated into the mod- 
els (e.g., Tozzi & Norman 2001; Voit et al. 2002; Oh & 
Benson 2003; Borgani et al. 2004; McCarthy et al. 2004, 
hereafter M04). 

Interestingly, by assembling published Chandra and 
XMM-Newton X-ray data from the literature, M04 un- 
covered a population of relaxed clusters (e.g., A2029, 
PKS0745, Hydra A) with central entropy profiles (where 

we define the "entropy" , S, as kTrie ^^'^) that are well 
approximated by powerlaws (see also the recent study 
of Piffaretti et al. 2004). M04 were able to account 
for these clusters (as well as other clusters that exhibit 
large entropy cores) with a simple analytic model that 
includes the detailed effects of cooling and that assumes 
the ICM was initially heated prior to cooling. In par- 
ticular, M04 found that if the ICM in a given cluster is 
initially heated^ by only a mild amount, cooling naturally 

^ We note that such "preheating" does not preclude subsequent 
heating events following cluster formation, such as AGN feedback 
initiated by the accretion of cold gas at the cluster center. 

^ However, the mild amount of heating is sufficient to delay 
the onset of catastrophic cooling for up to several Gyr and, as a 
result, the predicted global cold gas fractions remain consistent 
with observationally established limits. 



establishes approximate powerlaw entropy profiles near 
the cluster core^. It would be interesting to understand 
from a physical perspective what is driving this trend in 
the observed and model clusters. 

Unfortunately, a rigorous calculation of the effects of 
radiative cooling and the subsequent (quasi-hydrostatic) 
inflow of intracluster gas can only formally be obtained 
by numerically solving the time-dependent hydrody- 
namic equations (as approximately done, e.g., in M04). 
This makes it a challenge to understand physically why 
the above trends are established by cooling. Despite this 
inconvenience, there are, however, analytic tools at our 
disposal that can help to establish a physical picture. 
For example, Bertschinger (1989) (hereafter, B89) used 
a self-similarity analysis to derive the behavior of cooling 
flows in clusters. One of the interesting results derived 
by Bertschinger is that the logarithmic slopes of the gas 
density and temperature profiles in the limit r <g; rcooi 
(where rcooi is the radius at which the cooling time of 
the gas equals the age of the cluster) depend only on the 
shapes of the cluster gravitational potential and of the 
cooling function, A(r). The behavior of the gas density 
and temperature profiles can be used to infer how ra- 
diative cooling influences the distribution of intracluster 
entropy, which is perhaps a more fundamental quantity 
since convection will strive to prevent the establishment 
of a rising entropy profile towards the cluster center (see, 
e.g., Voit et al. 2002). Thus, the study of B89 is a good 
starting point for our investigation of the effects of cool- 
ing and inflow on intracluster entropy. 

In the present study, we briefly review the self-similar 
analysis of B89, including the basic assumptions made in 
that study and their validity, and use his results to de- 
rive how radiative cooling modifies the entropy profiles 
of clusters. We compare the self-similar solution to the 
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results of the 1-D coohng model of M04 for clusters with 
poworlaw dark matter profiles and that cool via thermal 
bremsstrahlung and, indeed, find extremely good agree- 
ment in the limit r <C Vcooi- We further demonstrate that 
the self-similar solution should provide an accurate de- 
scription for cooling (via both bremsstrahlung and line 
emission) in more realistic dark matter halos. Finally, 
we discuss a number of interesting uses of this result, in- 
cluding measuring the shapes of cluster gravitational po- 
tentials (that, e.g., could be compared against the usual 
hydrostatic equilibrium method), testing the reliability 
of numerical cooling routines in analytic models and hy- 
dro simulations, and as a simple way of setting up initial 
conditions for cluster models with radiative cooling. 

2. COOLING FLOWS AND SELF-SIMILARITY: THE 
BERTSCHINGER SOLUTION 

As noted above, in general, the time-dependent hydro- 
dynamic equations must be solved numerically in order to 
obtain a detailed picture of the effects of radiative cooling 
on cluster gas. However, when the time dependence is 
due to a single physical process that can be characterized 
by a unique scale length [in this case, the cooling radius^, 
fcoolit)], similarity solutions can provide useful physical 
insight. Adopting this approach, B89 derived the general 
behavior of cooling flows in clusters of varying gravita- 
tional potentials (and with varying cooling functions as 
well). We are especially interested in the properties of 
his solution in the limit r -C rcooi, since they are inde- 
pendent of the initial conditions of the gas. (Outside the 
cooling radius the initial conditions are of crucial impor- 
tance since radiative cooling hasn't had enough time to 
significantly modify the gas there.) Before examining the 
solution itself, let us first review the basic assumptions 
made in B89 and their validity. 

Since self-similar solutions can be characterized by only 
a single scale length, some simplifying assumptions are 
required in order to obtain a solution for the properties 
of cooling flows. In particular, Bertschinger assumed 
that the gravitational potential and the cooling func- 
tion could be characterized by simple powerlaws that 
remained fixed as a function of time. Of course, in real- 
ity, neither of these assumptions are strictly valid. High 
resolution numerical simulations indicate that the dark 
matter density profiles of clusters (and halos of other 
masses as well) have a characteristic scale length (the so- 
called scale radius, Vg), where the index of the powerlaw 
profile changes from relatively shallow (between ~ — 1 
and -1.5; e.g., Navarro, Frenk, & White (NFW) 1997; 
Moore et al. 1999) to relatively steep (^ —3). The cool- 
ing function is not scale-free either, owing primarily to 
line emission. Thus, one is justified in questioning the 
physical relevance of a model that invokes these assump- 
tions. (However, physical relevance may not be of major 
concern if one is simply testing the reliability of numerical 
cooling methods.) We examine in §3.2 whether relaxing 
these assumptions significantly affects the shapes of the 
resulting entropy profiles in the cores of massive clusters. 

^ The cooling radius grows larger with time. Unfortunately, 
the self-similar solution in B89 is expressed in terms of the initial 
cooling radius. Thus, it might be expected that comparison to 
observations, which are used to infer the present cooling radius, is 
somewhat ambiguous. However, as highlighted by B89, this leads 
to only a small error since, over the course of a cluster's life, the 
cooling radius grows only by a small amount. 



Other assumptions made in the analysis of B89 in- 
clude spherical symmetry, subsonic inflow, and single- 
phase cooling. Spherical symmetry is expected to be 
approximately valid, at least in an average sense for a 
reasonably large relaxed cluster sample. Likewise, sub- 
sonic flow should hold for the majority of the gas within 
rcooh except possibly near the very center where the gas 
may become transsonic. However, it has yet to be deter- 
mined whether or not multi-phase cooling is important in 
clusters. In the absence of significant non-gravitational 
heating, the intracluster gas is thermally unstable. How- 
ever, because the gas flows into the cluster center essen- 
tially as fast as it cools, we expect multi-phase cooling to 
be important only near the very center (see B89). Recent 
spatially-resolved Chandra and XMM-Newton spectra of 
the central regions cooling flow clusters have confirmed 
that, probably with the exception of the very central ra- 
dial bin, single-phase models provide at least as good a 
fit as multi-phase models (e.g., David et al. 2001; Mat- 
sushita et al. 2002). For the present study, we assume 
single-phase cooling. 

Finally, it is implicitly assumed that there are no sig- 
nificant sources of non-gravitational heating (such as 
AGN feedback, thermal conduction, and turbulent mix- 
ing) present in the ICM. In §4, we give a brief discussion 
of the potential impact of such heating. 

Implementing the above assumptions, B89 renormal- 
ized the hydrodynamic equations by removing any time 
dependence arising through rcooi{t)- A self-simflar so- 
lution is obtained if one neglects the acceleration terms 
(which is valid since the inflow of gas is highly subsonic) 
in the renormalized hydro equations. It is straightfor- 
ward to derive the limiting behavior of the gas density 
and temperature profiles in the limit r <^ Vcool under 
these conditions (see eqns. 2.30 of B89): 



dlogp _ -3-h(2-a)(l-/3) rflogT 



2- a (1) 



d log r 2 d log r 

where we have assumed that dark matter dominates the 
gravitational potential and that pdm and A(r) oc 

T^^. Note that the definitions of a and (3 differ from the 
definitions of these symbols in B89. 

In §1, we defined the "entropy", S, in terms of gas den- 
sity and temperature. The above equations can, there- 
fore, be used to yield the logarithmic slope of the entropy 
profile within rcooi- 
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(2 - a) + 1 (2) 



Thus, for (3 = 1/2 (i.e., cooling dominated by thermal 
bremsstrahlung) , 
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which yields 7 = 1 for a singular isothermal sphere {a = 

2), 7 w 1.4 for a Moore et al. profile in the limit r ^ 
Ts {a = 1.5), and 7 « 1.8 for a NFW profile in the 
limit T <C r.s {a = 1). Note, however, that the value of 
7 in equation (2) depends only weakly on the shape of 
the cooling function (that is, for reasonable values of /3 
ranging from —1/2 to 1/2) for a > 1. 

Below, we compare this simple analytic result with the 
1-D coohng model of M04. 
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3. COMPARISON WITH M04 
3.1. Powerlaw clusters 

M04 developed a simple radiative cooling code which, 
when applied to the entropy injection model clusters 
of Babul et al. (2002), successfully and simultaneously 
reproduces the luminosity-temperature and luminosity- 
mass relations (including their associated intrinsic scat- 
ter) and yields detailed fits to the entropy, surface bright- 
ness, and temperature profiles of clusters as inferred from 
recent high resolution X-ray observations. As alluded to 
in §1, these authors found that radiative cooling estab- 
lished a powerlaw entropy profile in cores of their model 
clusters that experienced only mild preheating (that is, 
for those clusters that were initially injected with < 300 
keV cm^, i.e., the entropy cooling threshold for massive 
clusters). It is interesting to see whether or not this trend 
can be explained by the self-similar solution of B89. 

In order to test this, we use the model of M04 (see 
§2.2 of M04 for a detailed discussion of the model) to 
track the effects of cooling for a set of clusters with ar- 
bitrary initial conditions (recall that the solution of B89 
within the cooling radius does not depend on initial gas 
conditions). For specificity, however, we show results for 
clusters that have a total dark matter mass of lO^^M©, a 
total gas mass of « 1 .5 x lO^^M© , and a maximum radius, 
fhaio, of 2.06 Mpc. The cluster gas is initially assumed to 
be in hydrostatic equilibrium within a dark matter halo 
that is characterized by a density profile pdm r^" (the 
normalization being set by the mass and radius given 
above). The cluster gas is then allowed to evolve via ra- 
diative cooling and inflow until a stable entropy profile 
is achieved. To calculate the cooling rate, we assume 
a cooling fimction that scales as A(T) oc T^/^ with a 
normalization that is set by matching a zero metallicity 
Raymond-Smith plasma at a temperature of T ss 10^ K. 
Finally, as in M04, we remove any gas that is able to 
cool below a threshold temperature of 10^ K from the 
calculation and assume that its dynamical effects on the 
cluster gravitational potential are negligible. Clearly, if 
a significant amount of gas is able to completely cool 
this assumption will be violated. However, in this case, 
we expect these effects will be relevant only for the cen- 
tral tens of kpe and will have only a minor effect on the 
overall entropy distribution within the cooling radius. 

To show how the effects of cooling are linked to the 
underlying gravitational potential, we present Figure 1. 
Focusing first on the top left-hand panel, we start with 
a cluster that is initially characterized by entropy profile 
that contains a core and a logarithmic slope of 7 = 1.1 
outside the core. This is the initial slope adopted by 
Babul et al. (2002) and M04 and is what one expects 
if the gas is in hydrostatic equilibrium and if its den- 
sity profile traces that of the dark matter (e.g., Voit et 
al. 2002; Williams et al. 2004), an expectation that is 
supported by high resolution hydrodynamic simulations 
(e.g., Lewis et al. 2000; Voit et al. 2003; Ascasibar et al. 
2003; Voit 2004). As can be clearly seen, the slope of the 
dark matter profile, a, is important in determining the 
resulting slope of the entropy profile, 7, within O.lrhaio 
(which corresponds roughly to Vcooi for these clusters). 
In particular, as the dark matter profile steepens the re- 
sulting entropy profile becomes more shallow, agreeing 
qualitatively and quantitatively with equation (3). 




0.01 0.1 1 0.1 1 

Fig. 1. — The effects of cooling and inflow on the entropy distri- 
bution of clusters. Solid lines represent the initial entropy profiles. 
The dotted, short dashed, and long dashed lines represent the re- 
sulting steady-state entropy profiles when the clusters are evolved 
in dark matter halos that have density profiles characterized by 
powerlaw indices of a = 2.0, 1.5, and 1.0, respectively. The var- 
ious panels show the resulting profiles for different initial entropy 
distributions. Figure 2 presents a comparison of the central loga- 
rithmic entropy slopes to the Bertschinger solution. 
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Fig. 2. — Comparison of the self-similar solution of B89 with 
the results of M04's cooling model. The shaded region reflects the 
uncertainty in the best fit powerlaw indices of the entropy profiles 
(within Tcool) shown in Fig. 1. 



The remaining panels of Figure 1, which show the re- 
sults for clusters that initially have steeper entropy pro- 
files outside the entropy core, illustrate that the initial 
gas conditions of the clusters have almost no effect on the 
resulting steady-state entropy profile within rcool- Like- 
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Fig. 3. — The effects of line emission and realistic dark matter profiles on the steady-state entropy profile. Left: A cluster with a Moore 
et al. dark matter halo. Right: A cluster with a NFW dark matter halo. In both panels the dotted line represents the initial entropy 
distribution, the thick solid lines represent the final entropy distribution assuming cooling with a pure thermal bremsstrahlung cooling 
function, and the thick dashed lines represent the final entropy distribution assuming cooling with a Raymond-Smith plasma model with 
Z = 0.3^0. The thin dashed lines indicate the slope predicted by B89's self-similar model assuming (3 = 1/2 and a = 1.5 (left panel) and 
a = 1.0 (right panel). They have been arbitrarily normalized to cross the thick lines at the cooling radius. 



wise, what happens in the interior of the cluster does not 

significantly influence gas outside of rcooi, as expected. 

We have fitted the entropy profiles within TcooI with 
powerlaws. However, we find that the entropy profiles 
within rcool are not exact powerlaws and there is some 
"wiggle" room in the best fit logarithmic slope, depend- 
ing on the range of radii over which the profiles are fitted. 
For example, the best fit logarithmic slope over the range 
< r < rcool /'^ differs slightly from the best fit over the 
range rcooi/'^ <r< rcool- We use these two radial inter- 
vals to roughly quantify the scatter in the best fit slope. 
Figure 2 presents a comparison between the analytic self- 
similar solution of B89 and the fits to the entropy profiles 
shown in Fig. 1. The shaded region roughly reflects the 
uncertainty in the best fit powerlaw indices for the pro- 
files predicted by M04's cooling model. 

Reassuringly, excellent agreement between the self- 
similar solution and the 1-D cooling code is obtained. 
Thus, the self-similar cooling wave model of B89 pro- 
vides a physical basis for the powerlaw trends found by 
M04. In addition, the agreement in Figure 2 gives us 
confidence in the reliability of the 1-D cooling model de- 
veloped in M04. 

3.2. Realistic clusters 

Observed clusters and clusters formed in cosmological 
numerical simulations do not have pure powerlaw gravi- 
tational potentials. Furthermore, the ICM contains a sig- 
nificant quantity of metals and, consequently, cools not 
only through thermal bremsstrahlung but also through 
line emission. Line emission has the effect of distorting 
the cooling function away from the powerlaw form that 
is characteristic of bremsstrahlung. For these two rea- 
sons, the physical relevance of the results presented in 
§3.1 may be questioned. Below, we investigate the ex- 



tent to which the shape of the resulting entropy profile 

is affected by these assumptions. 

Since we are introducing additional scales into the 
problem, it is important that we construct realistic clus- 
ter models. We consider two different systems: one with 
a NFW dark matter profile and one with a Moore et 
al. dark matter profile. Both systems have been cho- 
sen to have the same total gas and dark matter masses; 
specifically, Mgas{r2m) ~ l.SxlO^M© and Mdm{r2m) = 
IO^^Mq, where r2oo ~ 2.06 Mpc. We use a typical cluster 
dark matter concentration of Cnfw = ~ 3.4 for 

the NFW halo (e.g.. Eke et al. 1998; Bullock et al. 2001) 
and CMooro = r2m/r(-2) = cnfw/0.630 « 5.4 for the 
Moore et al. halo (see Keeton 2001). The above implies 
a scale radius, rg, of w 600 kpc. As for the intracluster 
gas, we turn to the studies of Voit et al. (2003) and Voit 
(2004). These authors foimd that the entropy profiles of 
a large sample of clusters generated with a numerical sim- 
ulation of a ACDM cosmology including hydrodynamics 
(e.g., shock heating) but not radiative cooling are ap- 
proximately self-similar over a wide range of masses (see 
Fig. 11 of Voit 2004). The initial entropy distributions of 
our model clusters are assumed to be identical to Voit's 
best fit to his simulated clusters^. The initial gas density 
and temperature distributions are determined through 
the equation of hydrostatic equilibrium by applying the 
boundary condition that the total amount of gas within 
r2oo is equal to that specified above. 

* At large radii, i.e., for r > 0.1r200, Voit (2004) reports a best fit 
entropy profile of S{r) oc r^'^. At small radii, however, there is an 
apparent entropy core whose origin remains uncertain. Wc have 
shrunk this core for computational convenience, since the model 
clusters reach steady state more quickly if they have small ini- 
tial entropy cores. However, this modification docs not affect our 
results or conclusions since the entropy core is contained entirely 
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To compute the eflFects of cooling, we again make use 
of the model developed by M04. In order to gauge the 
effects of line emission, we explore two different cooling 
functions: the pure thermal bremsstrahlung function im- 
plemented in §3.1 and a Raymond-Smith plasma with a 
metallicity set to 0.3Zq. As in the case of the powerlaw 
models, we neglect the dynamical effects of mass drop 
out and we run the cooling model until steady-state en- 
tropy profiles are achievc;cl. At steady state, we find that 
both clusters have similar global emission-weighted tem- 
peratures with kTew w 5 keV. 

In Figure 3, we plot the initial and final entropy pro- 
files of our model clusters. In both panels, the dotted 
lines represent the initial entropy distributions, the thick 
solid lines represent the final distributions when cooled 
using a pure bremsstrahlung cooling function, and the 
thick dashed lines represent the final distributions when 
cooled using the 0.3^0 Raymond-Smith plasma cooling 
function. The thin dashed lines indicate the slope pre- 
dicted by B89's self-similar model assuming (3=1/2 and 
a = 1.5 (left panel) and a = 1.0 (right panel). They 
have been arbitrarily normalized to cross the thick lines 
at the cooling radius. 

First, consider the role of the cooling function. For 
the Moore et al. halo, there is virtually no dependence 
on which cooling function we use. The resulting entropy 
distribution for the NFW halo, however, is slightly shal- 
lower if we include line emission. This difference can 
be understood as follows. The gravitational potential of 
the Moore et al. halo is steeper than that of the NFW 
halo. Consequently, gas flowing into the center of the 
Moore et al. halo requires more thermal support to re- 
main in hydrostatic equilibrium. In Figure 4, we plot the 
final (steady-state) temperature distributions of the two 
model clusters. Note that the central temperature of the 

within the cooling radius. As discussed above, the resulting steady- 
state entropy profile within Tcooi depends only on the shapes of the 
gravitational potential and the cooling function (and not on the 
initial properties of the gas within Tcooi)- 



Moore et al. halo is roughly 3 times larger than that of the 
NFW halo. In Figure 5, we show the cooling function for 
a Raymond-Smith plasma with Z = 0.3^0. For temper- 
atures of kT > 2 keV, the cooling function is dominated 
by thermal bremsstrahlung and is well approximated by 
a powerlaw; A(T) oc T^/^. This then explains why the 
entropy profile of the Moore et al. halo is unaffected by 
line emission: at any particular time there is virtually 
no gas below 2 keV (except at the exact center where 
gas rapidly cools below the X-ray emitting threshold of 
Rs 10^ K). The shallower NFW potential, however, per- 
mits some gas to cool below 2 keV before reaching the 
center of the cluster. The thin dotted line in Fig. 5 shows 
the best fit powerlaw to the cooling function between 0.1 
keV < kT < 2 keV, which has an index of /3 « -0.35. 
Using this value for /3 in equation (2), we are able to 
account for the slight 10%) deviation in the shape of 
the entropy profile within the central 0.01r2oo (~ 20 kpc) 
of the NFW halo. 

The shape of the final entropy profile is more sensitive 
to the shape of the gravitational potential than it is to 
the shape of the cooling function. Thus, we might expect 
self-similar solution to be a poor description of the final 
shape of entropy profiles in realistic dark matter halos. 
However, the left hand panel of Fig. 3 demonstrates that 
the shape of final entropy profile of the Moore et al. halo 
(thick lines) is virtually identical to that of a halo with 
a pure powerlaw profile of a = 1.5 (thin dashed line). 
Likewise, with the exception of the small deviation at 
the center due to line emission, the shape of the entropy 
profile in the NFW halo is essentially identical to that 
of a halo with a pure powerlaw profile of a = 1.0 (see 
right hand panel of Fig. 3). Recall that in the limit of 
r ^ Tg, the logarithmic slopes of the Moore et al. and 
NFW halos asymptote to a = 1.5 and 1.0, respectively. 
Fig. 3 illustrates that what is relevant is the shape of 
the local gravitational potential (i.e., at r < Tcooi)^ not 
the shape of the overall potential. Thus, the self-similar 
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kT (keV) 

Fig. 5. — The cooling function for a Raymond-Smith plasma 
with Z = 0.3Zq. The dotted line represents the best fit powerlaw 
(with /3 Ri -0.35) to the function over the range 0.1 keV <kT < 2.0 
keV. 



solution provides an excellent description of steady-state 
cooling entropy profiles in realistic clusters so long as the 
cooling radius is smaller than the cluster's scale radius. 
This condition should be met for most high mass clusters 
as the typical cooling radius of clusters is of order ~ 100— 
200 kpc (e.g., Peres et al. 1998), while the typical dark 
matter scale radius of massive clusters in high resolution 
cosmological simulations is of order ~ 400 — 700 kpc (e.g., 
Eke et al. 1998; Bullock et al. 2001). 

We conclude that the self-similar solution of B89 
should provide a good description of the shapes of en- 
tropy profiles (for r < rcooi) of massive clusters that 
cool via thermal bremsstrahlung and line emission and 
that have realistic dark matter profiles. The reason for 
this is that self-similarity is only mildly violated for high 
mass clusters. The deep gravitational potential wells 
of massive clusters ensure that most of the intracluster 
medium is quite hot {kT > 2 keV), even within the cool- 
ing radius, and thus thermal bremsstrahlung (which is 
scale-free) dominates the X-ray cmissivity. The poten- 
tial wells themselves have a characteristic scale radius 
but, for massive clusters, this radius is typically much 
larger than the cooling radius. Therefore, the central 
cooling flow essentially "feels" only a pure powerlaw po- 
tential (as assumed by B89). 

4. DISCUSSION 

Using the self-similar solution of B89, we have shown 
that radiative cooling and inflow lead to a characteristic 
entropy profile within the cooling radius that depends 
only on the shapes of the cooling function and the grav- 
itational potential. We have compared the self-similar 
solution to the cooling model of M04 and, reassuringly, 
find excellent agreement for clusters with powerlaw po- 
tentials and that cool via thermal bremsstrahlung. Fur- 
thermore, we have demonstrated that the self-similar so- 



lution is also valid for more realistic dark matter poten- 
tials and cooling functions that include line emission, at 
least for massive clusters where bremsstrahlung domi- 
nates line emission and the typical dark matter scale ra- 
dius is much larger than the cooling radius. This result 
explains the numerically-derived trends found in M04, 
which provide a good fit to a number of observed cooling 
flow clusters. Kaiser & Binney (2003) have also recently 
reported that cooling establishes a powerlaw trend be- 
tween ICM gas mass and entropy. This trend is also 
likely to be explained by the self-similar solution. 

The self-similar model of B89 implicitly assumes that 
there are no significant sources of non-gravitational heat- 
ing present in the ICM. Heating introduces an additional 
scale into the problem and, potentially, violates self- 
similarity. It is clear, however, that some form of heating 
has (or is) occurred in real clusters in order to prevent 
the so-called cooling crisis (see Balogh et al. 2001). Addi- 
tional evidence for heating comes from recent high reso- 
lution Chandra images of many cooling fiow clusters that 
reveal buoyantly-rising bubbles of hot plasma (e.g., Heinz 
et al. 2002; Blanton et al. 2003). These bubbles were pre- 
sumably blown by a central ACN and must be heating 
the ICM at some level. Other sources of heating, such 
as thermal conduction (e.g., Medvedev & Narayan 2001) 
and stirring due to the orbital motions of cluster galaxies 
(e.g., El-Zant et al. 2004). arc also a possibility. The rel- 
evant questions therefore are: (i) Is the level of heating 
sufficiently high to severely violate self-similarity? (ii) 
If so, is the heating distributed or restricted to only the 
very center of the cluster? (iii) In the event the heat- 
ing is episodic (such as AGN heating), when did the last 
heating episode occur (i.e., has enough time passed in 
order to approximately re-establish a self-similar cool- 
ing flow)? A number of relaxed clusters with published 
entropy profiles show evidence for large entropy cores 
(M04). Clearly, such systems cannot be explained by 
the present self-similar cooling model, as these clusters 
were severely heated and the heating was distributed 
to large radii. However, M04 (see also Piffaretti et al. 
2004) also found that several massive cooling flow clus- 
ters (e.g., A2029, PKS0745, Hydra A) have nearly pure 
powerlaw entropy profiles (except perhaps at very small 
radii, r < 30 kpc). This likely indicates that the self- 
similar model provides a reasonably accurate description 
of the cooling gas in this particular subset of observed 
clusters. 

A test of the above hypothesis is to infer the slopes 
of the gravitational potentials of these clusters by us- 
ing the slopes of their observed entropy profiles together 
with equation (3). This may then be compared with the 
results obtained using the standard hydrostatic equilib- 
rium method for inferring the gravitational potential of 
clusters. Making use of the clusters with powerlaw en- 
tropy profiles from M04 (A2029, PKS0745, Hydra A), 
we find that the logarithmic slopes of their total matter 
density (dark matter and baryons) profiles are relatively 
steep, with 1.3 < a < 2. This agrees quite well with the 
recent hydrostatic analysis of Chandra data of 10 relaxed 
cooling flow clusters by Arabadjis & Bautz (2004). Thus, 
for this small sample of clusters, the self-similar model 
appears to provide an apt description of the cooling gas 
in the centers of these clusters. A more detailed com- 
parison will soon be possible as the number of published 
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cluster entropy profiles is rapidly increasing. 

Quite independent of how well it describes observed 
clusters, the self-similar model also has a number of inter- 
esting theoretical uses. We briefly discuss but two here. 

1.) A simple method for calculating initial conditions 
of analytic model clusters with radiative cooling. This 
could serve as a "poor man's alternative" to a model 
that explicitly takes into account the effects of radiative 
cooling and inflow on intracluster gas. For example, a re- 
alistic set of initial conditions could be generated by using 
the results of non-radiative simulations (e.g., Lewis et al. 
2000; Loken et al. 2002; Voit 2004) to describe the gas 
at large radii (r > rcooi) while using the self-similar so- 
lution of B89 to describe the properties of the gas within 
Tcooi- The normalization of the entropy profile within 
fcooi (which is not specified by the self-similar solution) 
could be set by matching the non-radiative simulation 
results near Vcooi- One example of where such initial 
conditions might be useful is for models that explore the 
ability of various heating mechanisms to offset radiative 
losses of the ICM. For example, a number of recent AGN 
heating simulations (e.g., Quilis, Bower, & Balogh 2001; 
Ruszkowski & Begelman 2002) have invoked initial con- 
ditions such as isothermality or temperature profiles de- 
rived from non-radiative simulation (as opposed to that 
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expected for a cluster that is cooling radiatively) and this 
may have some effect on the estimates of the energetic 
requirements for the prevention of catastrophic cooling. 
It would be interesting to see whether the estimates of 
the amount of required heating change significantly when 
more realistic initial conditions (such as those proposed 
above) are implemented. 

2.) A test of the reliability of cooling routines im- 
plemented in analytic models and hydrodynamic simula- 
tions. Because the self-similar solution is a simple func- 
tion, it can easily be used to test, for example, how well 
various formulations of smoothed particle hydrodynam- 
ics (SPH) or mesh-based techniques [such as adaptive 
mesh refinement (AMR)] treat the effects of cooling (see, 
e.g., Abadi, Bower, & Navarro 2001). We are currently 
undertaking such a study using a variety of popular an- 
alytic and hydrodynamic codes (Dalla Vecchia et al. in 
preparation). 
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